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ABSTRACT 

We have calculated the emergent X-ray properties for models of cluster cooling flows including the 
effects of accumulated cooled material. The opacity of this cooled gas can reduce the overall X-ray 
luminosity of the cooling flow, and values of M based on these luminosities can underestimate the 
true value by factors of ~ 2. We find that accumulated cooled material can produce emergent surface 
brightness profiles much like those observed even for nearly homogeneous gas distributions. Consequently, 
much more of the gas may be cooling below X-ray emitting temperatures in the central regions of cooling 
fiows (r <^ 10 kpc) than one would infer from observed X-ray surface brightness profiles assuming the 
gas was optically thin. Similarly, the central densities and pressures in cooling flows may have been 
underestimated. We show that distributed absorption in cooling flows produces a number of observable 
effects in the spectrum which may allow it to be differentiated from absorption due to gas in our Galaxy. 
These include a characteristic suppression of the continuum below ~ 2 keV, absorption features such as a 
redshifted O K-edge, and diminished intensity of resonance emission lines. Spectra including the effects 
of intrinsic absorption are not well fit by foreground absorbing models. Attempting to fit such models to 
the spatially resolved spectra can lead to underestimates of the true absorbing column by factors of 3-20. 
Fits to integrated spectra of the entire cooling flow region can either underestimate or overestimate the 
mass of the absorbing gas depending on the specifics of the model. We discuss the potential detection 
of these effects with AXAF, XMM, and Astro-E. 

Subject headings: cooling flows — galaxies: clusters: general — galaxies: elliptical and lenticular, cD 
— intergalactic medium — radiative transfer — X-rays: galaxies 
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1. INTRODUCTION 

More than half of the clusters of galaxies observed at 
X-ray wavelengths exhibit evidence for cool gas in their 
cores ( Edge, Stewart, & Fabian 1992). The cooling time 
of this gas is often short compared to cluster ages; there- 
fore, the gas cools at rates that are often very large, 
Mc 10 — 2000 Mq yr~^. In the absence of a heating 
mechanism to balance cooling, gravity and pressure from 
the hotter, surrounding gas will drive the cooling material 
into the core of the cluster. If cooling flows are long-lived 
phenomena, these rates imply that ~ 10^^ Mq of material 
would cool over the lifetime of the cluster. Determining 
the flnal state of this cooling material remains the funda- 
mental problem concerning the cooling flow theory. 

The most obvious repositories for the cooling gas are 
cold molecular and atomic clouds, and stars. The cen- 
tral dominant galaxy in cooling flow clusters often have 
blue stellar colors in their central regions, which indicate 
that stars are currently being formed there (McNamara 
& O'Connell 1992). However, the observed star formation 
rates are generally 10% of Mc- Therefore star formation 
cannot account for the cooling material without appealing 
to unusual initial mass functions. Similarly, extensive ef- 



forts have been made to detect the accumulated cooled 
material either as gas at some temperature below the X- 
ray emitting regime (T <^ 10^ K). Gas is seen in cooling 
flows at lO'^ K (Heckman et al. 1989; Baum 1992) and in 
a few cases, as neutral atomic or molecular gas (Lazareff 
et al. 1989; Mirabel, Sanders, & Kazes 1989; McNamara, 
Bregman, & O'ConncU 1990; Jaffe 1992; O'Dea, Baum, & 
Gallimore 1994a). Dust is also seen in a number of cool- 
ing cores (Bregman, McNamara, & O'Connell 1990; Wise 
et al. 1993). In all cases, however, the detected levels of 
atomic and molecular gas are too low IO^-IO^^Mq) to 
account for the cooled gas which would accumulate over 
the age of the cluster. 

The detection by White et al. (1991) of excess X-ray 
absorption in a sample of cluster cooling flows was the 
first direct evidence for a sufficiently large mass of cold 
material. Using Einstein SSS (Solid State Spectrometer) 
spectra, these authors found that many cooling flow clus- 
ters exhibit significant levels of excess absorption over that 
expected from the Galaxy with typical excess column den- 
sities of 10^^ cm^^. Evidence for excess absorption in cool- 
ing flows has also been found in analyses of X-ray spectra 
taken with detectors on Einstein, ROSAT, BBXRT, and 
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ASCA (Lea, Mushotzky, & Holt 1982; Miyaji 1991; Allen 
et al. 1993; Fabian et al. 1994; Allen & Fabian 1997). The 
excess absorption columns detected by White et al. (1991) 
were only observed in clusters with spectral evidence for 
cooling flows and show some evidence for a correlation be- 
tween M and ANh- This correlation and the observed 
spatial coincidence between excess absorption and cool- 
ing flows suggests that the absorbing material is intrin- 
sic to the cluster and probably located within the cool- 
ing flow region. Integrating these excess column densities 
over the area of the cooling region in cluster cores implies 
the presence of large quantities of cold absorbing material 
(Mcoid ~ lO-'^-'^-lO^^M) and may represent the first direct 
evidence for the large amounts of cooled material which 
current coohng flow models predict (White et al. 1991). 

On the other hand, very extensive searches have been 
made to detect the excess absorber in emission or absorp- 
tion at radio wavelengths in lines of H I or CO and have 
not succeeded (e.g., McNamara & Jaffe 1993; Antonucci & 
Barvainis 1994; O'Dea et al. 1994b). It is becoming diffi- 
cult to understand how so much X-ray absorbing gas could 
have escaped detection in other wavebands (e.g., Voit & 
Donahue 1995). 

Often the ROSAT PSPC spectra of cooling flows are 
inconsistent with large columns of excess foreground ab- 
sorption (e.g., Sarazin, Wise, & Markevitch 1998), but are 
consistent with and may require large amounts of internal 
absorption in the cooling flow (Allen & Fabian 1997). For 
nearby clusters where the ROSAT or ASCA observations 
can resolve the cooling flow region, the excess absorption 
appears to be concentrated to the center of the cluster 
and cooling flow (r < 200 kpc) (Allen et al. 1993; Irwin & 
Sarazin 1995; Fabian et al. 1994; Allen & Fabian 1997). 

In the standard data analysis of X~ray spectra (e.g., in 
XSPEC or IRAF/PROS), it is conventional to treat ab- 
sorbers as lying in the foreground of emission components. 
This assumption allows the simple arithmetic combination 
of additive emission sources and multiplicative absorption 
components. However, X-ray observations suggest that 
the excess absorbing material in cluster cooling cores oc- 
cupies the same spatial region as the X-ray emitting gas 
(Allen et al. 1993; Irwin & Sarazin 1995; Allen & Fabian 
1997). Certainly, one would expect such a spatial cor- 
relation if the absorber originated as cooling X-ray gas. 
Including the effects of absorbing material which is inter- 
mixed with the X-ray emitting gas is not easily treated 
within the framework of conventional X-ray spectral mod- 
eling. Allen & Fabian (1997) used de-projection techniques 
based on ROSAT PSPC X-ray colors to attempt to deter- 
mine the three dimensional distribution of the absorption. 

In order to provide more accurate models of the X-ray 
spectra of cooling flows and to assess the effect of inter- 
mixed absorbing material, we have calculated the emer- 
gent X-ray properties for a set of inhomogeneous cool- 
ing flow models including the opacity due to accumulated 
cooled gas. For a given cooling flow model, we have self- 
consistently included the X-ray absorption due to vary- 
ing fractions of the total cooled material generated by the 
model. The details of the input models and the numerical 
solution of the radiative transfer equation are described 
in In we present X-ray luminosities, spectra, inte- 
grated surface brightness profiles, and emission line equiv- 



alent widths, and spectral profiles for the models. The im- 
plications of these results for the interpretation of existing 
X-ray observations and their application to observations 
with AXAF and Astro-E are discussed in §|[ In the final 
section (§||), we summarize our conclusions. We demon- 
strate that such large quantities of absorbing material can 
have a significant impact on the observed X-ray properties 
of cluster cooling flows. 

2. CALCULATIONS 
2.1. Inhomogeneous Cooling Flow Models 

As an input to the transfer calculations, we have se- 
lected a number of representative models for inhomoge- 
neous cooling flows, some of which were presented previ- 
ously in Wise & Sarazin (1993). These models were cal- 
culated using the numerical technique described in White 
& Sarazin (1987) and a subset have been used previously 
by Sarazin & Graney (1991) to estimate the optical coro- 
nal line emission from cooling flows. The optically thin 
X-ray properties for some of these models have been dis- 
cussed previously in Wise & Sarazin (1993, hereafter Pa- 
per I) while the effects of opacity due to cooling gas were 
discussed in Wise & Sarazin (1999, hereafter Paper II). 
We note that the opacity of gas which is in the process of 
cooling below X-ray emitting temperatures (as opposed to 
the accumulated cooled material described in this paper) 
was shown previously in Paper II to be insufficient to ac- 
count for the observed excess absorption. These cooling 
flow models have been discussed in detail in these papers 
and, in the interest of brevity, we will only briefly describe 
their relevant characteristics here. For a more detailed 
discussion, the reader is referred to these sources. 

The set of models discussed in Papers I and II span a 
range of gas deposition including both homogeneous mod- 
els as well as models with varying levels of inhomogeneity. 
Homogeneous models, by definition, exhibit a single den- 
sity and temperature at any given radius. It is important 
to realize that homogeneous cooling flow models cannot 
reproduce the observed properties of the excess absorp- 
tion. In a homogeneous model, all of the gas cools below 
X-ray emitting temperatures just within the sonic radius, 

, which is the point where the inflow velocity exceeds the 
sound speed in the gas (see §2.1 of Paper I and references 
therein). Values for the sonic radii are typically quite small 
(rg <^ 1 kpc); therefore, any accumulated cooled material 
would be confined to the very innermost regions of the 
cooling flow. In contrast, the spectral results of White et 
al. (1991) indicate that the absorbing material must cover 
a significant fraction of the cooling flow region (r ~ 100 
kpc), and the spatially resolved spectral measurements of 
Allen et al. (1993) and Irwin & Sarazin (1995) confirm 
this extended absorption. Consequently, in choosing mod- 
els for these calculations, we have only considered inho- 
mogeneous models which feature cooling material (and by 
implication cooled material) distributed over large spatial 
scales. 

The models used in this paper are listed in Table |l| and 
assume a spherically symmetric, steady-state inflow with 
only gravitational and gas pressure forces included. The 
outer boundary is taken to be the "cooling radius" , Tc, 
which corresponds to the point at which the instantaneous 
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isobaric cooling time, tc, is equal to the age of the clus- 
ter. In all cases, the cluster age was taken to be ta — 10^° 
years. As described in Paper I, the models are completely 
determined by three quantities: the total inflow rate M^; 
the gas temperature T^, and the mass deposition profile, 
p(r), which describes how inhomogeneous gas cools below 
X-ray emitting temperatures as a function of radius in 
the cluster. Both Mc and Tc are evaluated at the cooling 
radius Tc and values of Mc = 100 or 300 Mq yr~^ and 
Tc — 8.0 X lO'^ K were used for all models. 

To specify the mass deposition profile, p(r)^ two forms 
for the mass inflow rate were used. For most of the mod- 
els, the simple two-phase cooling flow model of White & 
Sarazin (1987) was employed. In this prescription, the 
mass deposition profile is parameterized as 



p(r) 



5'^ 



P{t) 



(1) 



where tc is the instantaneous isobaric cooling time, p{r) is 
the gas density, T{r) is the gas temperature, and P is the 
gas pressure. The total emission and cooling rate per unit 
volume in the gas is p^A{T), where A is the total emissiv- 
ity coefficient. The parameter q is referred to as the "gas 
loss efficiency parameter." The values of q for each of the 
models are given in column 5 of Table |l|. 

In addition to models using the White & Sarazin mass 
deposition formula (eq. [T]) , one model was included which 
used a form for the mass deposition profile derived by var- 
ious groups from deconvolutions of cluster X-ray surface 
brightness distributions. Fabian et al. (1984), Stewart et 
al. (1984), and Thomas et al. (1987) have argued on the 
basis of these observations that the inflow rate varies as 
M{r) oc r. In the C300_8_fb model, this assumption was 
adopted, which implies a mass deposition profile of the 
form 

Mr) ^ (2) 

We note that this model is quite similar to a White & 
Sarazin model with q w 3.4. 

The various physical parameters for each of the mod- 
els is given in Table 0. The models use a naming con- 
vention of Cxxxjy-zz, where xxx — Mc (in units of Mq 
yr~i), y = Tc (in units of 10^ K), and z = 10 x g. For 
the "Fabian et al." style model with M (x r, z — "fb." 
Three of the models were presented previously in Sarazin 
& Graney (1991) and Paper I, and the second column in 
Table |^ gives their identifications in the nomenclature of 
those papers. The remaining columns give the cooling rate 
Mc and temperature Tc at the cooling radius, the gas loss 
parameter, the cooling radius Tc, and the sonic radius Tj. 
A value of = indicates a fully subsonic solution. 

2.2. X~Ray Emission 

The calculation of the X-ray emission from the gas 
is discussed in considerable detail in Papers I and II. 
We included bound-bound, bound-free, free-free, and two- 
photon emission from ions of the elements H, He, C, N, 
O, Ne, Mg, Si, S, Ca, Fe, and Ni. Elemental abundances 
for these species were set at one half of the solar values as 
defined by Meyer (1979). The emissivities were calculated 



under the usual "coronal limit" conditions, as described in 
Paper I. Of course, the one major exception to these as- 
sumptions was the inclusion of the opacity of the gas. We 
assume that the cooled material is at temperatures be- 
low 10^ K, and produces no appreciable X-ray emission. 
At the column densities of interest, fluorescent line emis- 
sion by the cooled gas is not expected to be significant. 
However, if the columns of cold material in the absorbing 
clouds were larger, fluorescent emission of Fe K lines might 
be detectable (Churazov et al. 1998). 



2.3. X-Ray Opacity 

In calculating the opacity in the cooling flow, we have 
included contributions from the ambient hot gas, the cool- 
ing hot gas, and the accumulated cooled material. Of these 
three components, the first two (opacity due to ambient 
and cooling material) were discussed in detail in Paper II 
and are treated in the same manner here. Processes con- 
tributing to the opacity of the hot gas include resonant line 
absorption, photoelectric absorption, and electron scatter- 
ing. As discussed in Paper II, processes such as resonant 
scattering in the X-ray lines or electron scattering merely 
redistribute photons within the cluster and produce no net 
loss in the spatially integrated spectrum. We have also 
shown previously that the photoelectric opacity in an in- 
homogeneous cooling flow model due to ambient and cool- 
ing intracluster gas is small with r ^ 0.01 (see Paper II). 
Thus, the opacity from these gas phases cannot account 
for the observed excess soft X-ray absorption in cooling 
flows. 

In the present paper, we include a new source of opac- 
ity associated with the accumulated cooled gas. We as- 
sume that this material is relatively cold, so that its opac- 
ity is primarily due to photoabsorption. Values for the 
photoionization cross-sections were taken from an updated 
version of the Raymond & Smith (1972) emissivity code 
provided by John Raymond. The opacity of the cooled 
material will depend somewhat on its ionization state. We 
have assumed that the cooled material is cold and neutral. 
Although the physical state of the cooled material has not 
been established by observations. White et al. (1991), Fer- 
land, Fabian, & Johnstone (1993), and Daines, Fabian, & 
Thomas (1994) have argued that the absorbing material 
must be in the form of small, cold clouds to have avoided 
detection at other wavelengths. As mentioned above, the 
gas must be at a temperature < 10^'* for oxygen to retain 
K-shell electrons and thus present a significant photoelec- 
tric cross-section (Shull & Van Steenberg 1982). In fact, 
the total opacity should remain fairly constant for temper- 
atures much below this value, since the cross-sections are 
fairly insensitive to changes in the ionization state once the 
ion has regained a full complement of K-shell electrons. 

For a given model, we must also specify the amount 
and spatial distribution of the cooled material in order to 
calculate its contribution to the photoelectric opacity. As 
discussed in §|l|, our approach is to construct self-consistent 
models in which the absorbing material is assumed to rep- 
resent a fraction of all the gas which has cooled below X- 
ray emitting temperatures over the lifetime of the cooling 
fiow. Thus, we have taken the density of accumulated cold 
absorbing material as a function of radius in the cluster to 
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be 



Pabsir) = r] p{r) ta 



(3) 



where p{r) is the rate at which material is cooling below 
X-ray emitting temperatures and ta is the lifetime of the 
cooling flow, taken here to be 10^° years. The factor rj is 
a free parameter which allows us to vary the fraction of 
the cooling material which is stored as cold absorbing gas. 
Since it is essentially a deposition efficiency, we assume 
that 77 < 1. 

Inherent in the adoption of equation (|3|) is the assump- 
tion that the cold material accumulates where it cools. 
Obviously, this assumption is a simplification, since in 
practice a cloud of gas which is colder (and thus denser, 
assuming pressure equilibrium) than the ambient gas will 
tend to fall in towards the center of the cluster. White 
& Sarazin (1987) have calculated the distribution of the 
cooled material assuming it moves on ballistic orbits. Be- 
cause such material spends a significant fraction of its or- 
bital period at the largest distances from the cluster center, 
the resulting distributions are reasonably approximated by 
equation (O). Alternatively, Nulsen (1988) and Fabian et 
al. (1991) have argued that the cold clouds must comove 
with the background cooling flow. Such co-motion also 
results in a distribution similar to equation (^). In any 
event, the likely departures from equation (^) are in the 
sense that the profile of the cold gas is more centrally con- 
densed, which would tend to magnify the effects discussed 
in this paper. 

2.4. Numerical Solution of the Transfer Equation 

We have solved the transfer of radiation through the 
cooling flow numerically using the basic algorithm de- 
scribed by Yorke (1980; 1986; 1988). The transfer equa- 
tion is cast in the cylindrical (6, z) coordinate system of 
Hummer & Rybicki (1971). Note, we use b rather than 
p to represent the impact parameter which is fixed for a 
given "ray" through the flow and corresponds to a given 
projected radius on the sky. The coordinate z measures 
the depth into the cluster along each ray. With this choice 
of coordinates, the transfer equation can be written 



dlEjb, z) 
dz 



(4) 



where Ie is the intensity of radiation at photon energy E 
at any point in a given direction along a ray. The emis- 
sivity, eE^ contains contributions due to scattering and 
therefore depends implicitly on integrals of the intensity 
Ie- 

Having specified the variation of e^; and ke with radius 
within the fiow, we solve equation (|^) using the same nu- 
merical code employed in Paper II to calculate the effects 
of opacity due to the cooling gas. Spherical symmetry and 
a steady radiation field are assumed. The code is non- 
relativistic and assumes that velocity w, temperature T, 
and photon energy E satisfy \v\ <^ c, kT <^ nrieC^, and 
E <^ rrieC^. The effects of the fiow velocity due to the 
cooling flow have been included. This code is described in 
detail in Paper II and Wise (1992). For relevant details 
of the numerical technique, the reader is referred to these 
sources. 



3. MODEL RESULTS 

The models presented here are limited to the cooling 
flow region of the cluster (r < Tc). In practice, the cooling 
flow region is immersed within and smoothly connected 
to the general intracluster gas which extends out at least 
another order of magnitude in radius. We have chosen 
not to include the ambient cluster component to our mod- 
els in order to focus on the properties of the coohng fiow 
region and avoid introducing additional parameters asso- 
ciated with the cluster distribution. In practice, this de- 
cision means that some of the results of our models will 
differ from the equivalent properties directly observed in 
real clusters. For example, the X-ray surface brightnesses 
in our models drop to zero at the cooling radius (e.g., see 
Figure ^ below). Thus, to compare these profiles with ac- 
tual cluster surface brightness profiles, the emission due 
to the foreground and background cluster gas must be re- 
moved. 



3.1. Total X-Ray Luminosity 

Absorption due to the accumulated cold material will di- 
minish the total amount of X-ray emission from the cool- 
ing flow. The magnitude of this reduction in the total 
cooling fiow luminosity is shown in Figure |l| as a function 
of the deposition fraction rj. Curves are given for three 
cooling fiow models with a small degree of inhomogeneity 
q = 0.1 (C300_8_01), a moderate degree of inhomogeneity 
q — 1.0 (C300_8_10), and a high degree of inhomogeneity 
q — 4.0 (C300_8_40). Three different bandpasses are pre- 
sented, corresponding to the total X-ray flux (all photon 
energies > 0.1 keV), the flux observed with the ROSAT 
PSPC (Position Sensitive Proportional Counter) instru- 
ment, and the harder bandpass of the ASCA CIS (Gas 
Imaging Spectrometer). Also shown for comparison are 
the effects of a column Nh of foreground absorption. 

For values of the deposition fraction 77 ~ 1, the opacity 
due to the accumulated cold material can produce a signif- 
icant reduction in the total X-ray luminosity on the order 
of 40%, and an even larger reduction in the flux seen in 
the soft ROSAT VSPC bandpass (as large as 60%). The 
excess absorption columns seen by White et al. (1991) cor- 
respond to values of 77 w 0.1 to 0.4 (see Figure [l^ and §4.S 
below). We note that these smaller values of rj still result 
in appreciable reductions of '^30% in the total flux and 
-45% in the softer ROSAThand. 

The differences between the three models in Figure |^ 
are relatively small. At flrst sight, this lack of distinction 
is surprising, since a more homogeneous model {q = 0.1, 
C300_8_01) deposits more absorbing material in the center 
of the cluster where it covers less of the cooling flow re- 
gion. On the other hand, a more homogeneous model will 
also produce more of its X-ray emission in this region as 
well. In fact, equation (^ implies that there is always a 
correlation between the distribution of the X-ray emissiv- 
ity of the cooling gas (proportional to p) and the opacity 
of the absorption due to cooled gas. This correlation de- 
creases the model-to-model differences in the reduction in 
the X-ray luminosity. 



3.2. Emergent X~Ray Spectra 
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The primary evidence for substantial amounts of cold, 
absorbing material in the data of White et al. (1991) is a 
deficit of X-ray emission for photon energies below about 
2 keV over that expected from purely galactic absorption. 
To provide a basis for comparison with spectral observa- 
tions such as these, we have calculated the emergent X-ray 
spectra for our models. In Figure ||, the total integrated 
spectra from 0.1-10 keV are shown for a fairly homoge- 
neous model {q — 0.1, C300_8_01) for varying values of 
rj. An absorption component due to foreground Galactic 
material has not been included in this Figure. As Figure || 
demonstrates, the accumulated cooled material does pro- 
duce significant absorption below 2 keV even for moderate 
values of rj. 

For spatially resolved spectroscopy, the spectral signa- 
tures of excess absorption can be even more apparent. In- 
tegrated spectra average over the entire range of absorbing 
column densities in the cooling flow resulting in a lower 
mean level of absorption. However, a smaller aperture 
centered on the cluster will average over higher columns 
resulting in a greater mean absorption. Figure ^ shows sev- 
eral spectra for the same model extracted within a 10 kpc 
radius circular aperture centered on the core. The degree 
of absorption is correspondingly greater than in Figure |^. 

In addition to the level of soft X-ray absorption, it is 
also important to note that the accumulated cold mate- 
rial produces a distinctive signature in the shape of the 
emergent spectrum. Absorption due to foreground ma- 
terial produces an exponential fall off in the observed 
flux toward lower energies since, at a given photon en- 
ergy E, f{E) = fo{E)exp{—TE) where te is the optical 
depth. However, for absorbing material which is inter- 
mixed with the emitting gas, the dependence of the ab- 
sorption on optical depth is weaker. Roughly, we find that 
f{E) oc fo{E)/TE- Consequently for a given absorbing col- 
umn, foreground absorption will produce deeper absorp- 
tion edges and a stronger decrease in the spectrum toward 
lower photon energies than distributed absorption. 

This point is illustrated in Figure 0, which compares the 
spectrum for the g = 0.1 (C300_8_01) model with t] = 1.0 
(see Figure ||) with the spectrum of the same model with 
no internal absorption but foreground absorption corre- 
sponding to Nh = 2.5 X 10^^ cm~^. This amount of fore- 
ground absorption was selected to give roughly the same 
spectral shape as the internally absorbed model above 1.0 
keV. As noted, the model with foreground absorption de- 
creases much more steeply to lower photon energies, and 
has much deeper absorption edges (e.g., the O K-edge at 
0.532 keV). In principle, this difference in the shape of the 
emergent spectrum may provide a means to distinguish 
between foreground and intrinsic absorption. For clusters 
with large foreground Galactic absorption, howe ver , mak- 
ing such a distinction may be problematic (see §^^). 

3.3. X-Ray Surface Brightness Profiles 

Centrally peaked X-ray surface brightness profiles rep- 
resent one of the primary pieces of evidence for the ex- 
istence of cooling flows in clusters of galaxies. The X- 
ray emission due to the cooling gas is seen as an excess 
in the central surface brightness profile relative to that 
of the outer parts of the cluster (Jones & Forman 1984). 
In Figure |^, the emergent X-ray surface brightness pro- 



files are shown for a fairly homogeneous model (g — 0.1, 
C300_8_01) with values of r/ ranging from (no absorption) 
to 1.0 (100% of the accumulated cold material). These pro- 
files correspond to what would be observed in the ROSAT 
PSPC bandpass. Clearly, the central peak in the cooling 
flow X-ray emission is strongly depressed by the presence 
of distributed absorption. 

For comparison. Figure ^ also shows the X-ray surface 
brightness profile in the PSPC bandpass for the C300_8_fb 
model with M oc r and no internal absorption. This 
model provides a reasonable fit to the observed X-ray sur- 
face brightness profiles of many cooling flow clusters (e.g., 
Fabian et al. 1984). As Figure ^ illustrates, even if only 
10% of the cooling material collects in this cold form, sub- 
stantial reductions in Ix{b) will occur. The rough resem- 
blance of the intrinsically absorbed, essentially homoge- 
neous model profile to that of the optically thin, M (x r 
model highlights another important consequence of dis- 
tributed absorption. By suppressing the central peak of 
the surface brightness profile, a distributed absorber will 
make a given cluster cooling flow appear more inhomoge- 
neous than it actually is. 

ROSAT has, provided high quality X-ray surface bright- 
ness profiles for many clusters. However, much of the ini- 
tial work on the X-ray surface brightness profiles of cluster 
cooling flows was based on observations with the Einstein 
X-ray Observatory (e.g., Fabian et al. 1984). The Einstein 
Observatory had a substantially harder X-ray bandpass, 
extending to approximately 4 keV. Since the effects of ab- 
sorption (intrinsic or foreground) are lessened at higher 
energies, one might expect Einstein surface brightness pro- 
files to be less sensitive to the effects of distributed absorp- 
tion. To test this intuition and provide standards for com- 
parison to the Einstein data, we have also calculated the 
X-ray surface brightness for our models using the response 
of the Einstein IPC (Imaging Proportional Counter). Fig- 
ure ^ confirms that internal absorption has a smaller effect 
on Ix as seen by the Einstein IPC than ROSAT. Roughly 
speaking, one must increase rj by about a factor of two to 
produce the same effect on Einstein observations as with 
ROSAT. The qualitative effect on the surface brightness 
profile is the same. 

3.4. Resonant X-Ray Line Emission 

Strong X-ray line emission from cooling gas is one of the 
characteristic features of cluster cooling fiows ( Canizares, 
Markert, & Donahue 1988). Many of the X-ray lines 
produced by cooling plasma are non-resonance transitions 
having small oscillator strengths, /, and correspondingly 
small optical depths due to resonant absorption. Such 
lines are affected by the internal absorption in the same 
way as the adjoining continuum discussed above. Conse- 
quently, the equivalent widths of optically thin lines are 
not strongly affected by internal absorption. However, 
many of the strongest lines arise from resonance transi- 
tions and can be moderately optically thick in cooling 
flows (see Paper II). By itself, resonant line absorption 
(or scattering) merely redistributes photons spatially and 
spectrally. When coupled with a significant source of con- 
tinuum opacity, however, resonant scattering can increase 
the effects of the absorption on strong, resonance X-ray 
lines. This increase occurs for lines with significant optical 
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depth because the effective pathlength through the absorb- 
ing medium is increased. We consider some examples of 
these effects here. 

3.4.1. Line Equivalent Widths 

To illustrate some of the consequences of intrinsic ab- 
sorption for emission lines, we consider the observed prop- 
erties of the Fe XVII 2p^ - 2f3d resonance line 
which has a wavelength of 15.01 A {E ^ 0.826 keV). This 
line is produced by gas at temperatures in the range 1-3 
xlO^ K which is cooling rapidly out of the X-ray emitting 
regime. It has a large oscillator strength (/ = 2.35) and 
an optical depth due to resonant absorption of r ~0. 3-2.0 
depending on the assumed cooling flow model (see Pa- 
per II). Since this Fe XVII line has a low photon energy, 
it is also particularly susceptible to the effects of inter- 
nal soft X-ray absorption. In the simulations discussed 
here, we have assumed the line to be thermally broadened. 
This assumption maximizes the optical depth in the line. 
The gas in cluster cores, however, may be quite turbulent 
{vturb ~ 1000 km s^^) which would result in significantly 
broader emission line profiles and lower line optical depths 
(see Paper II). 

The reduction in the equivalent width, Wn. of the 
Fe XVII line as a function of rj is shown in Figure for sev- 
eral models with varying degrees of inhoniogeneity. Mod- 
els with q values of 0.1, 0.3, 1, and 4 are shown ranging 
from fairly homogeneous to strongly inhomogeneous, re- 
spectively. The model with Af cx r is also included for 
comparison. As Figure shows, resonant soft X-ray lines 
can be strongly attenuated by internal absorption in cool- 
ing flows. For a fixed value of 77, the equivalent width of 
the line decreases (relative to the optically thin value) as 
the deposition of the cooling material becomes more con- 
densed, i.e. as g decreases. Unlike the total luminosity of 
the flow, where the spatial coincidence of absorbing and 
emitting material serve to reduce model-to-model varia- 
tions, the increase in opacity due to resonant scattering 
is more sensitive to the location of the absorbing mate- 
rial. Thus, models with differing mass deposition profiles 
can produce very different reductions in the line fluxes. 
Note that the reduction in the equivalent width of the line 
in Figure is in addition to the overall reduction in the 
emission of continuum and optically thin lines at similar 
photon energies. Clearly, the fluxes of these lines can be 
greatly reduced. 

3.4.2. Line Surface Brightness Profiles 

Like their broadband counterparts, the presence of sig- 
nificant cold, absorbing material will suppress the central 
peaks in the surface brightness profiles of lower energy 
X-ray emission lines. Profiles for resonance lines will be 
particularly affected due to the interplay between absorp- 
tion from the accumulated cold material and the effects of 
resonant scattering. Figure ^ shows the emergent equiv- 
alent width of the resonance Fe XVII line as a function 
of projected radius for the fairly homogeneous q = 0.1 
(C300_8_01) model assuming various values of 77. As the 
amount of internal absorption increases, the gradient in 
the equivalent width flattens dramatically inside of r < 10 
kpc. Again, since the equivalent width compares the flux 



in the line to that of the adjoining continuum, the reduc- 
tion in equivalent width is in addition to the overall reduc- 
tion in the X-ray surface brightness shown in Figure ^ 

In the model without internal absorption (rj = 0), the 
surface brightness of the line is very strongly concentrated 
to the center of the cluster. This concentration simply 
reflects the fact that, for relatively homogeneous models, 
most of the gas cools inside the sonic radius which occurs 
at rather small radii (see Paper I). The Fe XVII line is 
emitted by cooling gas at relatively low temperatures, so 
the emissivity for this line is very large near the center of 
the model. Similar results are found for other resonance 
Fe L lines and for the resonance K-lines of low-Z elements. 
The K-lines of heavier elements (Si through Ni) are af- 
fected less by the opacity of accumulated cooled material 
due to their higher energy. 



3.4.3. Line Spectral Profiles 

The combination of continuum opacity from accumu- 
lated cooled material and resonance scattering will also al- 
ter the spectral line profiles of low-energy resonance lines. 
The left panel of Figure|| shows the spectral profile, (/>,,(-£), 
of the Fe XVII line in the integrated spectrum of a cooling 
flow model for different values of the efficiency parameter 
r/. These profiles were calculated for the fairly homoge- 
neous q = 0.1 model (C300_8_01) and assume thermal line 
broadening. The broad wings of the profile for the model 
without internal absorption are due to gas moving radi- 
ally inwards at a fairly high velocity near the center of the 
model (see Paper I). 

The right panel highlights the differences in the spectral 
profile due to optical depth effects and shows the ratio of 
the profiles with internal absorption to the profile with no 
accumulated cooled material. Overall, the combined ef- 
fects of resonance scattering and internal absorption serve 
to reduce the strength of the line. The broad wings of the 
line profile are produced by rapidly infiowing gas near the 
center of the model, where the opacity is largest. Conse- 
quently, the wings are reduced more than the line center, 
where there are contributions from larger radii. For lines 
produced by emitting material at larger radii or in models 
with smaller infiow velocities, it is the line center which 
is most strongly reduced. The spectral profiles depicted 
in Figure |^ also exhibit asymmetries relative to the line 
center. These asymmetries arise because blue-shifted pho- 
tons (due to material fiowing into the cluster center on the 
far side of the cluster) must traverse a greater absorbing 
column than red-shifted photons (due to material flowing 
into the cluster center on the near side of the cluster). 

Similar results are found for other soft X-ray resonance 
lines. Emission lines from higher ionization stages such 
as the Fe Ka at 6.7 keV are less effected by photoelectric 
absorption, although such lines can still be affected by res- 
onance scattering (Paper II). These effects are reduced if 
the lines are significantly turbulently broadened, as this 
reduces the line scattering optical depth. In any case, the 
detection of these effects will require very high spectral 
resolution observations; as a result, it may be some time 
before opacity effects on line shapes can actually be ob- 
served. 
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3.5. Absorption Edges 



In addition to its broadband effects (see §3^), accu- 
mulated cooled material in the cluster core will produce 
absorption edges in the emergent X-ray spectra of cooling 
flows. At the columns of interest here, the strongest spec- 
tral feature associated with the absorption is the O K-edge 
at 0.532 keV (e.g, Figure ^). If the intrinsic absorption is 
associated with the cooling flow, the O K-edge will be red- 
shifted to lower energies such that Eedge = 0.532/(1 -I- z) 
keV. Thus, two O K-edges will appear in the spectrum: 
an unredshifted edge due to Galactic gas (and possibly 
oxygen in the detector and spacecraft environment), and 
a redshifted edge. Detection of such a redshifted edge fea- 
ture at the redshift of the cluster would provide definitive 
evidence that the excess absorbing material was associ- 
ated with the cooling flow and not due to anomalously 
high foreground absorption. Arnaud & Mushotzky (1998) 
have reported detection of such a redshifted O K-edge in 
the Perseus cluster based on analysis of BBXRT spectra. 

We can measure the observed depth of the O K-edge 
from our self-absorbed cooling flow models directly. In 
Figure |lO|, the optical depth in the edge is plotted for a 
range of deposition efficiencies, rj, in the various models. 
The observed optical depth was determined assuming the 
decrement in the observed fiux at the edge energy was 
given by exp(— Tedge)- Note that this assumption is equiv- 
alent to assuming a foreground absorption model. For 
?7 = 1, the measured depth in the edge ranges between 
r ~ 0.1-0.3 for all the models. Given an absorption cross 
section at O K of tr ~ 6.0 x 10^^^ cm^, these optical depths 
imply columns of ANh ~ 2-5 x 10^° cm^^ assuming the 
absorbing materia l is in the foreground of the cluster. As 
we show below in §4.3, this "foreground" assumption leads 
to significant underestimates of the true excess absorbing 
column. No foreground Galactic column was assumed in 
Figure |l^. 

The spatial variation of the O K-edge optical depth is 
shown in Figure For each model, the cooling region was 
divided into annular bins 5 kpc wide and the optical depth 
in the edge was then determined for the spectrum from 
each bin. The different spatial profiles apparent in Fig- 
ure ^ simply trace the degree of inhomogeneity, and hence 
the degree to which the absorbing material is centrally con- 
densed, in the underlying model. For a given model, Tedge 
increases steadily as one moves to smaller projected radii 
until the point where the cooling flow becomes optically 
thick in the edge. In the model with AI{r) oc r, this point 
occurs at a radius of ~ 50 kpc; while in the slightly in- 
homogeneous g = 0.1 model, Tedge increases steadily into 
a radius of ~ 10 kpc. Since this "saturation" point cor- 
responds to the radius at which the actual optical depth 
in the edge is r 1, Figure |ll| implies that the measured 
value of Tedge Can be used to estimate the true absorb- 
ing column at this radius. From Figure |l|, the observed 
values for Tedge at these saturation radii range between 
~ 0.25-0.30 implying underestimates of the true columns 
by factors of ^ 3-4. This values are consistent with the 
results presented below based on spectral fits to the in- 
tegrated spectra assuming foreground absorption models 
(see §U). 

In practice, it is difficult to separate an intrinsic, red- 
shifted O K-edge from a Galactic edge with the ASCA 



SIS (Sohd Imaging Spectrometer) or AXAF ACIS {AXAF 
CCD Imaging Spectrometer) because of their spectral res- 
olution at these energies. However, such observations may 
be possible with AXAF for some higher redshift clusters 
(z ^ 0.1), where the intrinsic edge is separated from the 
Galactic feature by ^ 50 eV. Although the High Energy 
Transmission Grating (HETG) spectrometer on AXAF 
has high spectral resolution, its throughput is low and its 
resolution is degraded when observing extended sources. 
Astro-E, on the other hand, should easily be capable of 
detecting redshifted absorption edge features in cooling 
fiow spectra. 

3.6. Scaling of Model Results 

Obviously, it is not possible to present results cover- 
ing the full range of all cluster cooling flow parameters in 
this study. Thus, we have given a few characteristic ex- 
amples to demonstrate the effects of internal absorption. 
However, it is generally possible to crudely assess the ef- 
fects of internal absorption for other parameter values by 
approximate scaling of the models presented here. The 
steady-state models we have presented are primarily char- 
acterized by flve parameters: the total cooling rate Mc, 
the temperature at the cooling radius T^, the age of the 
cooling flow ta] the gas loss efficiency parameter q; and 
the fraction of accumulated cooled gas which is present in 
an X-ray absorbing form rj. In addition, the models de- 
pend on the parameters characterizing the gravitational 
field of the cluster and central galaxy, and the abundances 
in the cooling gas. However, if the temperature at the 
cooling radius is sufficiently high and one is not concerned 
with the details of the innermost parts of the flow near the 
sonic radius, the hydrodynamical solutions do not depend 
strongly on the detailed form of the gravitational poten- 
tial. The atomic physics of the emission processes does 
depend in a complex fashion on the abundances and the 
outer temperature, so we will assume these are fixed in the 
scaling given below. 

For a fixed temperature Tc and fixed abundances, the 
density pc at the cooling radius Tc depends inversely on 
age, pc (X t~^. The total cooling rate in a model scales 
approximately as Mc oc pcr^/ta- Thus, the cooling radius 
in a model with fixed abundances and Tc varies approxi- 
mately as 

Tc^M'jHl/'. (5) 

Equation ^ implies that the total mass of the cold ab- 
sorber in the model is 



Mats = V ta . 



(6) 



Consequently, the column density of the cold absorber 
Nniabs) and the implied absorption optical depth Tabs 
both scale as 

Tats (X Nniabs) oc ^ ocvM^H-'/' . (7) 

Intervening absorption affects the observed properties in 
nonlinear way; therefore, the value of Tabs must be the 
same for two models if they are to be scaled to one an- 
other. 

On the other hand, most of the emission properties of 
cooling flow models nearly scale with the total cooling rate. 
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Thus, the luminosity in any spectral feature (and the flux 
in that feature in the integrated spectrum) approximately 
scale as 

oc . (8) 

The surface brightness of most features will scale approx- 
imately as 



L 



(9) 



Resonance emissi on l ines present an exception to these 
scaling relations (§ |3.4[) . The properties of these lines de- 
pend on both the optical depth due to absorption by cold 
material Tabs and due to resonance scattering Tj.g. The 
latter scales approximately as 



CX PcTc OC Ml/3 . 



(10) 



Requiring this quantity to be constant would generally 
make it difficult to scale the models. In what follows, we 
will consider the scaling of all properties of the models ex- 
cept those associated with resonance emission lines which 
are optically thick to resonance scattering. 

Let us assume one wishes to determine the properties 
of a model characterized by Tc, q, Mc, ta, and rj which is 
not presented in this paper. Further, assume that a model 
with the desired Tc and q is presented, but with a differ- 
ent total cooling rate M° and age Then, from the grid 
of models presented in this paper, one should select the 
model with a value of the deposition efflciency parameter, 
?7°, given by 




-1/3 



(11) 



We note that this scaling can break down for models with 
large values of 77 and significantly different cooling rates or 
ages. Strictly speaking, equation |ll| could produce values 
of 77° which are much greater than unity. This much ab- 
sorption could not be produced by material cooling at the 
present rate. 

Once the value of if has been fixed, the luminosity and 
integrated spectrum of the desired model can be scaled 
from the corresponding fiducial model using 



while the surface brightnesses scale as 



(12) 



T K T° 




1/3 



-4/3 



(13) 



We have verified these relations for a number of cooling 
flow models, and flnd that they work to within 10-20% in 
essentially all cases. 



4. EFFECTS OF INTRINSIC ABSORPTION ON DERIVED 
QUANTITIES 

With the exception of Allen & Fabian (1997), most pre- 
vious analyses of the X-ray spectra from cluster cooling 



flows have not considered the effects of additional absorp- 
tion due to accumulated cooled material. Traditionally, 
any absorption in the cooling flow spectra has been as- 
sumed to come from foreground material and in many 
cases has been flxed at the Galactic value as measured 
by Stark et al. (1992). Neglecting the effects of this ma- 
terial, if present, can lead to serious errors in the derived 
physical properties of the gas in the cluster core. In the 
following discussion, we examine the effects of this inter- 
mixed internal absorption on a number of physical quan- 
tities commonly derived from X-ray observations. 

4.1. Total Cooling Rates 

Assuming the X-ray luminosity from a cooling flow is 
due to material cooling isobarically from some tempera- 
ture T, to flrst order Lx is simply proportional to the 
enthalpy of the cooling gas: 



5 kT 
2 firUp 



(14) 



For a measured total luminosity, Lx, this expression can 
be inverted to obtain an estimate of the total cooling rate 
of material. As we demonstrated in §3.1, even moder- 



ate amounts of accumulated cooled material can produce 
significant reductions in the observed X-ray luminosity. 
Therefore, estimates of the total cooling rate Mc based 
on the observed value of Lx will tend to underestimate 
the true value if the effects of internal absorption are ne- 
glected. Based on the results in Figure ^ the derived val- 
ues of Mc might be a s much as a factor of ~2 too low. As 
we show below (§4.3), estimates of the X-ray cooling rate 
based on fits to X-ray spectral observations which assume 
only foreground absorption will tend to underestimate the 
true absorbing column and therefore Mc as well. 

Individual X-ray emission lines can also be used to mea- 
sure the total cooling rate of material in cluster cores. If 
the line is emitted primarily from material at tempera- 
tures below the ambient cluster temperature Tc, then the 
luminosity in the line is given by 



5,v k 

-Mc 

2 /imp 



A;.„e(r) 
A(T) 



dT 



(15) 



where AundT) is the emissivity coefficient of the line. 
For lines emitted by gas at temperatures which are smaller 
than Tc, the luminosity is nearly independent of Tc, the de- 
tails of the cooling, or even the abundances (which affect 
both Aiine{T) and A(T) in a similar manner). This tech- 
nique requires high spectral resolution observations and 
has been applied to a number of clusters based on FPCS 
data from the Einstein X-ray Observatory (Canizares, 
Markert, & Donahue 1988). The High Energy Transmis- 
sion Grating (HETG) onboard AXAF should also provide 
data from several cluster cooling cores with sufficient spec- 
tral resolution to apply this technique. 

For optically thin X-ray emission lines, the effects of in- 
trinsic absorption will reduce the observed luminosity in 
the line, Lune, to the same degree as the nearby contin- 
uum. Thus, lines at lower energies will suffer greater atten- 
uation and consequently underestimate the total cooling 
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rate, Mc, to a greater degree relative to lines at higher en- 
ergies. For example, assuming a value of 77 ~ 1, optically 
thin lines with energies near 0.6 keV could suffer reduc- 
tions as large as 60% whereas the luminosity of a line at 2.0 
keV would be reduced by only 20%. These reductions are 
in addition to any reduction due to the intervening Galac- 
tic absorbing column. The actual observed reductions will 
of course depend on the true underlying deposition profile 
of cooling material and the detector response. If the line 
is a resonance line with an appreciable optical depth, an 
addit iona l reduction of about a factor of two is possible 
(see § 3^ and Figure 0) . 

As an aside, we note that differing amounts of intrinsic 
absorption in cluster cooling cores may be partially re- 
sponsible for the scatter observed in correlations between 
Mc and other cooling flow properties. Given the many po- 
tential sources of scatter in observed cluster correlations, 
it seems unlikely that accumulated cold material could be 
solely responsible. However, as these models illustrate, the 
magnitude of these effects is such that the presence of this 
material could represent a significant contribution to the 
observed scatter in such correlations. 



4.2. Spatial Distribution of Cooling Material 

Centrally peaked X-ray surface brightness profiles are 
one of the characteristic features of cluster cooling flows. 
However, it has been known for some time that the ob- 
served profiles are less centrally peaked than would be ex- 
pected if the gas in cluster cooling flows was homogeneous 
and all of the gas flowed into the inner regions (r ^ 1 
kpc) before cooling (Thomas et al. 1987; White & Sarazin 
1987). The conventional explanation of this discrepancy is 
that the gas is very inhomogeneous, and that denser lumps 
of gas cool more quickly. Thus, a significant portion of 
the gas cools below X-ray emitting temperatures before it 
fiows in very far, and does not contribute to the central 
X-ray emission. One way of quantifying this model is to 
calculate the amount of hot gas M(r) which is still flow- 
ing inward at each radius r. By mass conservation, this 
quantity is equal to the amount of gas which cools within 
this radius in steady-state. Thomas et al. (1987) have ar- 
gued that many clusters are reasonably well represented 
by a mass deposition proflle of the form M[r) oc r out to 
r ^ 200 kpc. One of our models (C300_8_fb) incorporates 
this spatial variation in the cooling rate. 

The implied form of the mass deposition profile M{r) 
will depend on the observed X-ray surface sur face bright- 
ness profile Ix{b)- However, as discussed in f|3.3|, intrin- 



solid line. Clearly, a value of 77 « 10% is sufficient to bring 
the surface brightness of this nearly homogeneous g = 0.1 
model into rough agreement with the optically thin M oc r 



model. In §4.3 below, we show that this level of internal 



sic absorption from accumulated cooled material has the 
effect of reducing the central surface brightness of cool- 
ing flows. Consequently, the gas in cooling flows may be 
less inhomogeneous than a simple interpretation of the ob- 
served surface brightness profiles might imply. This ef- 
fect is illustrated in Figure || which shows the predicted 
ROSAT PSPC X-ray surface brightness profiles for a 
slightly inhomogeneous q — 0.1 model with different values 
of 77 (expressed as percentages). As mentioned previously, 
cooling flow models with M oc r produce surface bright- 
ness profiles in reasonable agreement with X-ray observa- 
tions. For comparison, the optically thin surface bright- 
ness profile for the M oc r model is indicated by the heavy 



absorption would correspond to an observed excess column 
of about 1 X 10^^ cm~^ if the spectrum were fit assuming 
only foreground absorption. This value is similar to the 
excess columns which have been derived from X-ray spec- 
tra of cooling flows under this assumption (White et al. 
1991). 

As noted above, the variation of the inflow or cooling 
rate M(r) as a function of radius has been used to charac- 
terize the degree of inhomogeneity of a cooling flow. We 
have treated the predicted ROSAT X-ray surface bright- 
ness data from our models as observed data, and depro- 
jected the emission to determine M{r). Figure [ij shows 
the deprojected mass deposition profiles corresponding to 
the ROSATX-iay surface brightness profiles shown in Fig- 
ure ^. For comparison, the heavy solid line shows the ex- 
pected profile for a model with M(r) cx r (C300_8_fb). As 
the degree of internal absorption is increased (77 ~> 1), the 
derived M{r) profiles steepen which would be interpreted 
as evidence that more of the gas is cooling below X-ray 
emitting temperatures at larger radii. Again, a value of 
77 « 0.1 would cause the surface brightness of this fairly 
homogeneous g = 0.1 model to have a mass deposition 
profile which resembled M(r) oc r. 

These qualitative results apply to all of the models. We 
have illustrated the case for the fairly homogeneous g = 0.1 
model because the effects are particularly striking here. 
However, we note that the specific values of rj required to 
turn a more centrally condensed mass deposition profile (or 
equivalently, surface brightness profile) into one which re- 
sembles the observed M(r) oc r will depend on the specifics 
of the input model and the detector response. For exam- 
ple, due to its harder bandpass, if the above comparison 
were made with Einstein IPC data, an rj value a factor of 
two larger would be required to produce the same effect. 
Similarly, a much smaller spectrally determined excess col- 
umn of order 3 x 10^" cm~^ would convert a moderately 
inhomogeneous q — I model into a M(r) oc r model. 

As a final point, we note that the above effects do not 
apply to completely homogeneous models. In such mod- 
els, the gas does not cool to low temperatures until it is 
inside the sonic radius which typically occurs at very small 
radii (r 1 kpc). Therefore, any intrinsic absorption due 
to this accumulated cooled material would also be concen- 
trated within a very small region. If present, this material 
would produce a large reduction in the surface brightness 
profile but only at the very center of the cluster. It would 
not affect the surface brightness at larger radii since no 
absorbing material would have accumulated there. 



4.3. Absorbing Columns Estimated from Spectral Fits 

The evidence for excess absorption in cluster cooling 
flows has come mainly from spectral fits such as those by 
White et al. (1991). As we demonstrated in fi3.2|, absorb- 



ing material intermixed with the emitting plasma produces 
a characteristic spectral shape in the emergent spectrum 
which is substantially different from the simple exponential 
fall-off created by foreground absorption. Consequently, 
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fits to spectra from self-absorbed cooling flows assuming 
only foreground absorption will yield erroneous values for 
the derived columns. Specifically, such fits will tend to 
underestimate the true absorbing column. To quantify 
this effect, we have fit foreground absorption models to 
the emergent X-ray spectra of our models with internal 
absorption. In the following analysis, we have used the 
effective area and spectral resolution of the ROSATPSPC 
due to its good soft X-ray response. 



4.3.1. Fits to Integrated Spectra 

Most observations of excess absorption in cluster cool- 
ing flows have been based on the integrated spectra for the 
entire cooling flow region. For example, the original study 
by White et al. (1991) utilized the Einstein SSS detec- 
tor with a fleld of view of 6 arcmin. Most nearby cluster 
cooling flows are several arcmins in radius which means 
that the SSS gives the integrated spectrum of the entire 
cooling flow. Similarly, ASCA has a resolution of several 
arcmins and typically yields integrated spectra for the en- 
tire cooling flow region. The ROSAT PSPC has a better 
angular resolution of about 0.5 arcmin, but still only gives 
a few separate spectra within the cooling flow regions of 
nearby clusters (e.g., Allen et al. 1993). For comparison 
to such integrated spectra, we have determined the best- 
fit foreground column AN (spectral) using the equivalent 
integrated spectra from our models. 

To perform the spectral fits, we have used the following 
procedure. For a given cooling flow model and choice of 
deposition efficiency, 77, we have calculated the integrated, 
emergent spectrum for the entire cooling flow. This spec- 
trum was then convolved with the instrument response of 
the ROSAT PSPC to obtain the number of "observed" 
counts per energy bin, F"^" for this model and choice of 
rj. The instrument response consists of the effective area, 
Aeff{E), and the redistribution function, R{E,PHA), 
produced by the PSPC's finite spectral resolution. As a 
fitting model, we used the optically thin (77 — 0) spectrum 
for the selected cooling flow multiplied by varying amounts 
of foreground absorption and convolved with the PSPC in- 
strument response. For a given choice of Nh, this model 
yields the expected counts per energy bin, f/^*{Nh)- The 
value of Nh which yielded the best fit was then obtained 
by minimizing the statistic 



E 



[Fr-FriNn)] 



(16) 



where i is an index assigned to the spectral energy bins. 
The J^/** in the denominator of equation 16 assures the 
correct behavior if the errors in the observed spectra are 
dominated by Poisson errors from photon noise. Finally, 
because all clusters are observed through some amount of 
Galactic material, both and f/'* were also multi- 

plied by varying amounts of Galactic foreground absorp- 
tion prior to fitting. Fits assuming Galactic columns of 
0.1, 0.2, 1.0, and 2.0 x 10^^ cm~^ were calculated and these 
columns were subtracted from the best-fit foreground ab- 
sorption values after the fit. 

Figure |l^ gives the excess column densities determined 
using the above procedure as a function of the deposition 



efficiency rj. Using this Figure, we can compare the proper- 
ties of our models with a given 77 to the observed properties 
of clusters for which the excess absorption ANh {spectral) 
has been derived. Note that the values in Figure |l^ 
are specific to Mc = 300 Mq yr"^ and ta = 10^° yr. 
The scaling to other values is given in equation 0. The 
typical excess columns densities found by White et al. 
(1991) assuming foreground absorption models are of or- 
der ANh (spectral) ~ 1 — 2 x 10'^^ cm^^. These column 
values are in the regime where all of the model curves in 
Figure ^ cross, and correspond to « 0.1 to 0.4 for all of 
the models. 

Observationally, values of ANh (spectral) have been 
used to derive limits on the total amount of cold, X- 
ray absorbing material in cooling flows. Consequently, it 
is useful to compare the spectrally derived columns with 
those actually present. The sense of the difference de- 
pends on whether one considers spatially resolved spec- 
tra or the integrated spectrum of the entire cooling flow 
region. For spatially resolved spectra, the spectrally de- 
rived column assuming only foreground absorption is al- 
ways smaller than the actual column along that line of 
sight. This underestimate occurs because, for intermixed 
absorption, a fraction of the absorbing material always 
lies behind a portion of the emitting gas, and thus is not 
fully effective in producing absorption. However, if one 
considers the integrated spectrum of an extended region 
where the absorbing column and emission vary with pro- 
jected position, the difference can have either sign. A given 
amount of absorber is more effective if it is concentrated 
where the emission is most intense, and such a concentra- 
tion may not be resolved in the integrated spectrum. 

The spectrally-derived excess columns in Figure |l^ are 
determined from the integrated spectrum of the entire 
cooling flow region. Following equation 0, the average col- 
umn density of hydrogen in our models over the cooling 
flow region is given by 



ANh (true) 



lA7m„TTr^ 



lAlmpirr^ 



(17) 



Figure |lj compares the actual average column densi- 
ties determined using equation O with those derived from 
fits to the integrated spectrum of the cooling flow region 
(Fig. O). The heavy solid line indicates ANH(spectral) = 
ANh (true). For low column in the inhomogeneous 
models, the spectrally-derived columns overestimate the 
real column. For large columns of intrinsic material, 
ANH(true) > 1 x 10^^ cm~^, and for the most inhomo- 
geneous models, the spectrally derived column is smaller 
than the real average column for all of the models by a 
factor <^ 4. For the more inhomogeneous models (q > 2) 
the measured columns underestimate the true values down 
to columns of ANh (true) ~ 1-2 x 10^*^ cm^^. 

These results can be compared to a similar analysis per- 
formed by Allen & Fabian (1997) for a set of cooling flow 
clusters observed with the ROSAT PSPC. Using XSPEC 
simulations, these authors have simulated the errors in the 
apparent excess column for two distributions of internal 
absorption: a partial covering model and a "multilayer" 
model which assumes a homogeneous distribution of ab- 
sorbing throughout the X-ray emitting region. These in- 
put models were then fit with a simple uniform foreground 
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absorption model similar to our procedure. Of the two, 
this latter model is the closest approximation to the de- 
tailed distributions of accumulated cooled material mod- 
eled here. For values of ANnitrue) of 1 and 5 x 10^^ 
cm~^ and a Galactic column of Nnigal) of 1 x 10^° cm~^, 
these authors find that the fitted column densities under- 
estimate the true "multilayer" columns by factors ~ 4-5. 
These results are in good agreement with the factors of 
~ 3-4 obtained from our models. Allen & Fabian (1997) 
also note that the degree of underestimation is essentially 
independent of the Galactic column, a result we confirm. 

For the more homogeneous models 9 < 2, the spectrally- 
derived column can actually over-estimate the true average 
column for small amounts of absorption. This behavior is 
seen in Figure |l^ for values of AA^/f (true) ^ 5 x 10^° 
cm~^. In these models, both the emission and absorp- 
tion are strongly concentrated to the center of the cooling 
flow. When the emitting and absorbing material are con- 
centrated in this way, a smaller amount of absorber is more 
effective. One way of understanding this effect is to note 
that in fairly homogeneous models, the effective radii con- 
taining most of the emission and absorption are both much 
smaller than the cooling radius. Thus, if one averages the 
absorption over the projected area of the entire cooling 
flow, as in equation |l7|, the resulting value of the aver- 
age column density ANnitrue) will be much smaller than 
the true value in the region where the absorbing material 
resides. 



4.3.2. Fits to Spatially Resolved Spectra 

With the exception of a few nearby clusters observed 
with ROSAT, the cooling cores of clusters have not been 
well resolved spatially by instruments with good spec- 
tral resolution. Consequently, studies of excess absorp- 
tion have depended primarily on analyses of integrated 
spectra. However, AXAFs excellent spatial resolution 
(FWHM ~ 0.5 arcsec) will provide spatially resolved spec- 
tra for many cooling flows. With such future data in 
mind, we have calculated the values of AN (true) and 
ANH{spectral) as a function of projected radius in the 
cooling flow. To determine AN h (spectral), we divided 
the cooling region into annular bins 5 kpc wide. This 
choice gives roughly 20 bins over the cooling radius for 
each model. For a choice of rj, the emergent spectra within 
annular bins was calculated from the model surface bright- 
ness using 



~dE 



{bi:b2) = 2TT db' b' lE{b') (18) 

Jbi 



where bi and 62 represent the inner and outer radii of 
the bins. Similarly, emergent spectra were extracted from 
the corresponding optically thin model for the same set of 
annuli. These spectra were then fit using the same pro- 
cedure outlined above to calculate the best-fit values of 
ANh (spectral) . 

For comparison, the true variation of ANnib) was calcu- 
lated directly from the distribution of accumulated cooled 



material, Pabs(f) (see equation ||), using 
1 /■ 2 



ANH(b) 



lAlfJlr, 



Pabs dl 



1.47m 



""dr 'P^''^''^ 



P Jb 



(19) 

where / is the distance along the line of sight at a pro- 
jected radius b. Given AiV//(&), the average column den- 
sity within an annular bin was then taken to be 



27r / db' b' ANH(b') 

< ANh > (h : 62) = -^'^ . (20) 

7r(62 - 61) 

We note that equation ^ reduces to equation ^ in the 
limit that 61 — » and 62 — > re- These expressions were 
used to derive the true average column densities for the 
same annular grid used in determining AN [{(spectral). 

Figure 15 shows the results of this spatial analysis for 
two models: the slightly inhomogeneous q = 0.1 model 
(C300_8_01) and the model with M(r) oc r (C300.8_fb). 
As expected, the spectrally derived column density un- 
derestimates the true column at all projected radii. The 
magnitude of this error increases as one moves towards the 
center of the cluster. For the model with M(r) cx r, the 
column determined from spectral fits varies over a range 
of ~ 3-10 as one moves to smaller radii. In the g = 0.1 
model, more of the absorbing material is concentrated to- 
ward the center of the flow producing underestimates of 
the true absorbing column by factors of ~ 2-30. As before, 
the Galactic column has been subtracted from the curves 
in Figure nsl AXAPs spatial resolution should yield spa- 
tially resolved spectra, and consequently measurements of 
ANh, on spatial scales comparable to this example for 
clusters out to redshifts of z ^ 0.1. 



5. CONCLUSIONS 

In this paper, we have calculated models for the X-ray 
properties of cooling flow clusters including the effects of 
opacity due to accumulated cooled material. For consis- 
tency, we have assumed that the internal absorber has the 
same distribution as the gas which is cooling below X-ray 
emitting temperatures in the models. We have character- 
ized the amount of accumulated absorbing material as a 
fraction 77 of the amount of gas which would cool over the 
hfetime of the cluster (ta = 10^° yr) at the present rates 
of cooling. The emergent properties of the X-ray emis- 
sion were then calculated including the effects of radiative 
transfer through this accumulated material. 

The presence of accumulated cooled material will reduce 
the observed X-ray fluxes and the implied cooling rates of 
cooling flow clusters below the actual values. Depending 
on the cfhciency with which cooling gas is converted into 
absorbing material, the fluxes and cooling rates may be un- 
derestimated by as much as a factor of two over the 0.1-10 
keV energy range. For instruments with softer bandpasses, 
such as ROSAT, the reduction can be higher. Although 
modeling the effects of the intrinsic absorption with fore- 
ground absorption models will reduce the magnitude of 
this underestimate, these models tend to underestimate 
the total amount of absorbing material and therefore do 
not completely compensate for this effect. 
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Intrinsic absorption produces a characteristic spectral 
morphology below ^ 2 keV which is not well fit by fore- 
ground absorption models. Specifically, internal absorp- 
tion does not increase in effectiveness as an exponential 
function of the optical depth, as foreground absorption 
does. As a result, the integrated X-ray spectra for the 
models described here do not decline as rapidly at softer 
energies and the absorption edges are less deep than for 
models with foreground absorption. For spatially resolved 
spectra toward the core of the cooling flow, this effect is 
enhanced. We find that these differences should be dis- 
cernible in present ASCA and future Astro-E and AXAF 
spectra. On the other hand, it is difficult to distinguish 
spectra of models with internal absorption from those with 
foregroimd absorption with a partial covering factor. 

In addition to its effect on the broadband X-ray spec- 
trum, accumulated cooled material will produce redshifted 
absorption edges such as the O K-edge at an energy of 
0.532/(1-1-2:) keV. These redshifted edges provide a direct 
means of determining whether the absorbing material is 
intrinsic to the cooling flow and not due to anomalously 
high foreground absorption. Resolving such redshifted O 
K-edge features from edges due to foreground material and 
O in the CCD detectors of ASCA and AXAF is diflicult 
for nearby clusters given the spectral resolution of these in- 
struments. These features may be detectable in more dis- 
tant clusters with AXAF and XMM where the redshifted 
O K edge is more cleanly separated or in nearby cooling 
flows with Astro-E which has greater spectral resolution. 

The effects of internal absorption can have important 
implications for the analysis of X-ray surface brightness 
profiles from cluster cooling flows. Because the absorbing 
material in our models is concentrated toward the central 
regions of the cooling flow, it acts particularly to dimin- 
ish the central peak in the surface brightness and flatten 
the observed proflles. Starting with Einstein observations, 
it has been noted that cluster cooling flows are less cen- 
trally peaked than would be expected if the cooling gas 
were homogeneous. This fact has led to the canonical view 
that the gas in cooling flows is rather inhomogeneous and 
that much of the gas cools at large radii. However, our 
models indicate that opacity due to accumulated cooled 
material can produce emergent surface brightness profiles 
much like those observed even for nearly homogeneous gas 
distributions {q = 0.1). Consequently, much more of the 
gas may be cooling below X-ray emitting temperature in 
the central regions of cooling flows (r <J 10 kpc) than one 
would infer from observed X-ray surface brightness pro- 
files assuming the gas was optically thin. Another way 
of presenting the same result is to note that internal ab- 
sorption steepens the mass deposition profile. M{r), one 
would derive from the observed surface brightness distri- 
bution, making a nearly homogeneous model resemble one 
with M(r) oc r. We note that the reduction in the cen- 
tral surface brightness due to internal absorption would 



also mean that central cooling times have been overesti- 
mated, and that central gas densities and pressures have 
been underestimated. 

Low energy resonance X-ray emission lines are par- 
ticularly sensitive to the presence of internal absorption. 
These lines can have appreciable resonance scattering op- 
tical depths in cooling flows, and this scattering enhances 
the chances of absorption by accumulated cooled material. 
For lines such as the Fe XVII line at 0.826 keV, the inter- 
play between continuum opacity and resonant scattering 
can lead to reductions of 20%-80% in the equivalent width 
of the line. Like their broadband counterparts, estimates 
of the cooling rate based on such lines will underestimate 
the true cooling rate. Line surface brightness proflles and 
spectral profiles arc likewise effected. With sufficient spec- 
tra resolution, one can resolve distortions in the spectral 
profiles of these low energy resonance lines due to absorp- 
tion. However, detection of such signatures requires spec- 
tral resolution well beyond that currently available. Our 
models assume thermal line broadening. If significant tur- 
bulence is present in the cooling flow gas, the magnitude 
of these effects is reduced. Emission lines from high ion- 
ization stages, such as the Fe Ka line at 6.7 kcV, arc rela- 
tively unaffected by internal absorption though they may 
still have appreciable resonance scattering optical depths. 

Previous studies of intrinsic absorption in cooling flows 
have primarily determined the amount of absorbing ma- 
terial by fitting foreground absorption models to the ob- 
served X-ray spectra. As our models indicate, the spec- 
tral signatures of accumulated cooled material are not well 
fit by such models. For spatially resolved spectra, the 
spectrally derived columns always underestimate the ac- 
tual columns along the line of sight, typically by a factor 
of ~ 3-20 depending on the specifics of the cooling flow 
model. If the amount of excess absorption is derived from 
integrated X-ray spectra for the entire cooling flow (as has 
generally been the case to date due to the limited spatial 
resolution of the spectrometers used), the situation is more 
complicated. For large columns or fairly inhomogeneous 
models, the absorption is always underestimated this way. 
However, for small columns in nearly homogeneous mod- 
els, the degree of excess absorption can be overestimated. 
This effect occurs because both the emission and absorp- 
tion in nearly homogeneous models is strongly centrally 
concentrated, and a small amount of absorbing material is 
more effective in this case. 
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Table 1: 



Cluster Cooling Flow Models 



Model 


Previous 


Mc 


Tc 


Gas Loss 


T 






Model 


{Mq yr-i) 


(K) 




(kpc) 


(pc) 


C300.8_01 




300 


8.0 X 10' 


q = 0.1 


95.6 


582 


C300.8_03 




300 


8.0 X 10^ 


q = 0.3 


93.5 


414 


C300.8_10 


C3 


300 


8.0 X 10^ 


q = 1.0 


87.6 


115 


C300.8_40 


C4 


300 


8.0 X 10^ 


Q = 4.0 


72.3 





C300.8Jb 


C5 


300 


8.0 X 10^ 


M (xr 


83.7 


11 


C100.8_10 




100 


8.0 X 10^ 


q = 1.0 


61.3 


69 
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Fig. 1. — The percentage reduction in the X-ray luminosity as a function of the efficiency parameter r]. The top panel 
corresponds to the total luminosity from 0.1-10 keV. The middle and bottom panels show the reduction in the ROSAT 
PSPC and ASCA GIS bandpasses, respectively. In the three panels, curves are shown for several different inhomogeneous 
models considered in the text. The solid line corresponds to a model with a moderate q value of 1.0 while the dashed 
and dotted lines correspond to the extremes of the model grid. The dashed curve shows the reduction for a model with 
q = OA while the dotted line shows a model with q = 4.0. The heavy solid line in each panel indicates the percentage 
reduction in the various luminosities for a simple foreground absorption model with varying column densities Nu- The 
abscissa for the foreground absorption curves is indicated along the top axis of the plot. 
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Fig. 2. — Total integrated X ray spectra as a function of photon energy for the q = 0.1 (C300_8_01) model. Curves are 
shown for differing values of the efficiency parameter, rj, described in the text. A value of ry = corresponds to the the 
optically thin case. Each curve is labeled with the corresponding value of rj. 
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Fig. 3. — The central X-ray spectra as a function of photon energy for the q — 0.1 (C300_8_01) models. This is the 
spectrum from a central circular aperture with a projected radius of 10 kpc. The notation is the same as Figure ^ 
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Fig. 4. — The emergent X-ray speetrum from the q = 0.1 (C300_8_01) model with intrinsic absorption given by ry = 1 is 
compared to the same cooling flow model with no internal absorption, but with a foreground absorption oi = 2.5 x 10^""^ 
cm~^. Note that these spectra have similar shapes near 1 keV, but the foreground absorption model falls off more rapidly 
at lower photon energies. 
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Fig. 5. — The emergent X-ray surface brightness profile in the ROS AT PSPC bandpass as a function of projected radius 
for different values of the efficiency parameter, ij. Curves are shown for the nearly homogeneous g = 0.1 model and rj = 
corresponds to the optically thin surface brightness profile. The heavy solid line represents the expected surface brightness 
profile for a model with M(< r) oc r. Each curve is labeled with the corresponding deposition efficiency rj expressed as a 
percentage. 
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Fig. 6. — The emergent X-ray surface brightness profile as measured in the harder Einstein IPC bandpass. The notation 
is the same as Figure 1^. 
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Fig. 7. — The percentage reduction in the total equivalent width, Wr,, of the resonance Fe XVII line at 0.826 keV as 
a function of the deposition efficiency r] discussed in the text. Reductions are calculated relative to the optically thin 
equivalent width, Wo, for that model. Curves are shown for varying values of g: q — 0.1 (C300_8_01, dash-dot line), 
q = 0.3 (C300.8_03, dash-dot-dot-dot line), q = 1.0 (C300_8.10, solid hne), and q = 4.0 (C300_8_40, dotted line). The 
reduction for a model with M oc r is also included for comparison (C300_8_fb, dashed line). 
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Fig. 8. — The equivalent width, W^, of the resonance Fe XVII line at 0.826 keV as a function of projected radius for the 
fairly homogeneous q = 0.1 (C300_8_01) model assuming various values of rj. Individual curves are marked with the value 
of T] assumed. The thick solid line shows the profile for the M ocr model with no internal absorption for comparison. 
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Fig. 9. — The left panel shows the integrated spectral line profile, (pniE), for the Fe XVII line at 0.826 keV for different 

values of the efficiency parameter rj. These profiles correspond to a fairly homogeneous q = 0.1 (C300_8_01) model, and 
assume thermal line broadening. The continuum flux has been subtracted from each profile. The right panel shows the 
ratio of the emergent profiles to the model with no internal absorption. Curves arc shown for values of t] = 0.0 (solid 
line), T] = 0.01 (dotted line), r] = 0.1 (dashed line), and r] = 1.0 (dash-dot-dot-dot line). 
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Fig. 10. — The optical depth in the O K edge at 0.532 keV in our models is shown for different values of the deposition 
efficiency parameter rj. No foreground Galactic absorption has been included. Curves are shown for models with q = 0.1 
(dash-dot line), q = 0.3 (dash-dot-dot-dot line), q = 1.0 (solid line), q = 4.0 (dotted line), and M ocr (dashed line). 
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Fig. 11. — The optical depth in the O K edge at 0.532 keV is shown as a function of projected radius for our models. 
The model spectra were accumulated in 5 kpc annular bins and a value of 77 = 1 was assumed. The notation is the same 
as in Figure l9. 
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Fig. 12. — The deprojected mass deposition profiles, M{r), derived from the ROSAT PSPC X-ray surface brightness 
profiles calculated for the q = 0.1 model (C300_8_01) and shown in Figure]^. Each curve is labeled with the corresponding 
value of 7] expressed as a percentage. For comparison, the heavy solid line shows the expected mass deposition profile for 
a model with M(< r) oc r. 
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Fig. 13. — The excess absorbing column density for the various models derived from spectral fitting assuming foreground 
absorption is shown as a function of "q. The ordinate values represent the measured excess column after having been 
folded through the ROSATVSPG response. For all models, a foreground absorbing column of Nh = 2.0 x 10^° cm~^ was 
assumed and has been subtracted. All of the models have = 300 yr~^ and Tc = 8 x 10^ K. Curves are shown 
for models with q = 0.1 (dash-dot line), q = 0.3 (dash-dot-dot-dot line), q = 1.0 (solid line), q = 4.0 (dotted line), and 
M cx r (dashed line). 
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ANnltrue) [10'' cm"'] 

Fig. 14. — A comparison between the actual average internal absorption columns and the values measured from fits to the 
spectra. The actual average excess column is shown on the abscissa and is given by ANnitrue) — Mctai] / {1 Almpirr'^) . 
The ordinate values are the same spectrally derived columns shown in Figure O. The models and notation are the 
same as in that Figure, except that the heavy solid line indicates a one-to-one correspondence between ANnitrue) and 
AN (spectral). 
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Fig. 15. — The variation in the true average internal absorption column and the fitted column as a function of projected 
radius for two cooling flow models. The slightly inhomogeneous q = 0.1 model (C300_8_01) model is denoted by square 
symbols while the model with M{r) oc r (C300_8_fb) is denoted by the stars. The model spectra were accumulated in 5 
kpc annular bins and fit with a foreground absorption model as in Figures |l^ and ^ A value of 77 = 1 was used in both 
cases. For a given model, the solid line shows the true average column in the bin and the dashed line shows the resulting 
fitted values. A foreground absorbing column of Nh = 2.0 x lO^*' cm^^ was assumed and has been subtracted. 



